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Abstract — The objective of this work is, in the first place, the 
integration in a fusion process using hybrid DSmT model, both, 
the contextual information obtained from a supervised ICM clas- 
sification with constraints and the temporal information with the 
use of two images taken at two different dates. Secondly, we have 
proposed a new decision rule based on the DSmP transformation 
to overcome the inherent limitations of the decision rules thus 
use the maximum of generalized belief functions.The approach 
is evaluated on two LANDSAT ETM+ images, the results are 
promising. 

I. Introduction 

The management and the follow-up of the rural areas 
evolution are one of the major concerns for country planning. 
The satellite images offer a rapid and economic access to 
accurate homogeneous and updated information of studied 
territories. An example of application which results from this 
is related to the topic of the changes cartography, in this paper, 
we are interested to study the most subtle changes of the Argan 
land cover and other themes in the region of Agadir (Morocco) 
by contextual fusion /classification multidates based on hybrid 
DSmT model [1]— [3] and ICM with constraints [4], [5]. 

Our work environment, is the theory of Dezert-Smarandache 
[1]— [3] which is recent and very little implemented or used 
before the covered work of this paper, it was applied in mul- 
tidate fusion for the short-term prediction of the winter land 
cover [6]— [10] and, recently, for the fusion and the multidate 
classification [1 1]— [14], although the theory of evidence, it is 
more exploited for fusion/classification [11], [12], [15]— [19] 
also, for classifier fusion [20]-[22]. 

Our methodology can be summarized as following, after 
preprocessing of the images, a supervised ICM classification 
with constraints [4], [5] is applied to the two images, in order 
to recover the probabilities matrices for an after using in a step 
of fusion/classification basing on the theory of plausible and 
paradoxical reasoning known as Dezert-Smarandache Theory 
(DSmT) which allows to better assign the suitable pixels to 
the appropriate classes and also to detecte the changes. 

In this paper, in section 2, we describe in further detail 
the mathematical basis of the recent theory of plausible and 



paradoxical reasoning (DSmT), and give a description of our 
decision rule. The background and the details of the ICM 
algorithm with constraints are presented in section 3. In 
section 4, we provide the results of our experimentation where 
the algorithm was applied to a LANDSAT ETM+ image, the 
classification results are discussed in the same section followed 
by conclusions in section 5. 

II. Dezert Smarandache theory (DSmT) 

A. Principles of the DSmT 

The DSm theory was conceived jointly by Jean Dezert and 
Florentin Smarandache [1]— [3], it is a new way of representing 
and fusioning uncertain information. DSmT, considered as 
a generalization of the evidence theory of Demspter- Shafer 
[15], was developed to overcome the inherent limitations 
of DST (Dempster-Shafer Theory ) [l]-[3], [6]-[9], [13], 
[15]. The basic idea of DSmT rests on the definition of 
the hyper power set , from which the mass functions, the 
combination rules and the generalized belief functions are 
built. 

The hyper power set D G is defined as the set of all 
composite propositions/subsets built from elements of 0 with 
U and D operators such as: 

We define hyper power set , D G as follow : 

1) <p,9 1 ,02, ...On e D e 

2) If X,Y e D e , then A U Y £ D e and A n Y G D e . 

3) No other elements belong to D e , except those obtained 
by using rules 1) and 2) [l]-[3]. 

with 0 = ...6n},(/) is empty set. 

We define a map as follows: 

m s (.) :D 0 ^ [0,1] (1) 

associated to a given body of evidence s as 

mM)= 0 (2) 

and 

Y m s( x ) = 1 ( 3 ) 

xeD e 
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with m s (X) is called the basic belief assignment/mass (bba) 
of X made by the source s. 

The DSmT contains two models : the free model and the 
hybrid model [?], [1], [2], [7], [9], [11], the first presents limits 
concerning the size of the hyper power set D & , whereas 
the second has the advantage of minimizing this size, for this 
reason, it will be used in the continuation of our study. 

B. Combination rules 

The masses from the sources must be combined using a 
combination rule to have a new masses distribution for the 
elements of the hyper power set in order to promote an item 
compared to others. 

Within the framework of DSmT, there are several rules 
combination [1]— [3], [13], in our application, we have applied 
and implemented the majority of these rules in order to choose 
those which allow us to have good performances such as the 
PCR5 (conflicts proportionally with respect). 

C. Generalized belief functions 

From the function of basic mass, the generalized belief 
functions are defined , which model the imprecision and the 
uncertainty according to the hypothesis considered by a given 
source. 

The generalized belief functions used in this study namely 
the Credibility (Cr), Plausibility (PI) and DSmP Transfor- 
mation are defined for D 0 in [0, 1] and are given respectively 
by: Belief, plausibility and DSmP of an element X E D & \ 



Bel(X) = ^ 

vCX 

V e D e 

pim = £ m(v) 

vnx^f 

veD @ 



DSmP e (q b) = 0 
VX E G @ {f} , DSmP e (X) = 



< 



E 



V m(Z)+e-C(XnY) 

^ zcxnY 

C(Z )= 1 

E ZQy rn(Z^C(Y) 

C(Z) = 1 



m(Y ) 



( 4 ) 

( 5 ) 



( 6 ) 



where e > 0 is an adjustment parameter and D e possibly 
reduced by the introduction of the integrity constraints of 
the hybrid DSmT model. C(X fl Y) and C(Y) respectively 
indicate the cardinalities of the Y D X and Y. 



D. Proposed decision rule 

The decisions after combinations could be taken from 
generalized belief functions, to decide the belonging of a pixel 
to a given class, two cases are distinguished: 

• The pixel belonging to a simple class (used to improve a 
classification): in this case we use one of the following 



decision criteria: maximum of bba, maximum of the 
Credibility ( Cr ) (with or without rejection), maximum 
of the Plausibility (PI), Appriou criterium, (DSmP 
criterium • • • etc ) [6], [7], [11], [13] . 

• The pixel belonging to a composed class(e.g. in the case 
of change detection), in this case we cannot use the 
functions quoted previously because they are increasing 
functions and unsuited to the decision for the elements 
of union and of intersection. 

In an original step, we have proposed a new decision rule 
based on DSmP transformation and confidence interval to take 
in account the composed classes. In this decision rule we 
exploited the confidence interval [Bel (X )) , PI (X)\ by the 
definition of a new measurement which we have named : 
Global uncertainty IncG which is the sum of the uncertainties 
(Inc) of the hyper power set elements: 

VX € D e Inc (X) = PI (X) - Bel (X) (7) 

This new decision rule described in Algorithm 1 is applied as 
follows: 

For a given pixel x, we compare the Global uncertainty 
IncG of this pixel with a threshold (chosen experimentally). 
If it is lower than this threshold, the pixel is affected to the 
simple class which maximizes the DSmP transformation of all 
the simple classes, if not, it is affected to the composed class 
which maximizes the mass (gbba) of all the composed classes. 



Algorithm 1 The proposed decision rule 

l: IncG = E xeD & (PK X ) ~ Bel(x)) 

2: if (IncG < threshold) then 
3: if DSmP [x\ (0i) = max {DSmP [x\ (Si) 

with 1 < i < n then 
4: x E 0i 

5: end if 

6: else 

7: if m [x] (nOi) = max {[x] (fl Of) 

with 1 < i < \D°\ — (n + 1 ) then 
8: X E nOi 

9: end if 

10: end if 



III. ICM CLASSIFICATION WITH CONSTRAINTS 

The information contained in a satellite image is usually 
in the form of homogeneous objects. Indeed, an image of 
rural areas often consists of large homogeneous parcels, and 
therefore, an acceptable classified image must respect this 
property. Thus, the use of Markov Random Fields (MRF) takes 
in account this property of the neighborhood influence of a 
pixel on it. and therefore insists on coherence between the 
class of a pixel and that of its neighbors. It is a powerful 
mathematical tool for regularizing the classification of the 
satellite images. 

Moreover, the Markovien formalism constitutes a gateway 
to introduce a several constraints (spatial context, map of 




contours, temporal context, etc), for this reason, we have used 
the suggested method by [4], [5] as a method of classification 
in order to generate the probabilities for DSmT. This technique 
of classification provides, in addition to the constraint of 
regularization, a new constraint of segmentation so to refine 
the classification, these contextual constraints are controlled 
by a parameter of temperature in an iterative algorithm of 
optimization ICM (Iterated Conditional Mode). 



B. Trainings samples 

The samples are created automatically by using the 
ENVM.O software. We have identified 5 topics of land oc- 
cupation: Argan ( A ), Built /Oued (BO), Vegetation ( V ), 
Greenhouses ( Gh ) and Bareground (Bg) while basing on 
images already classified and on the ground truths, the figure 
2 represents their spectral responses.. 



IV. Result and discussion 
A. Study area and used data 

The study area is located in the region of Souss, it is in 
southern Morocco. Geologically, it is the alluvial basin of 
the Oued Souss, separated from the Sahara by the Anti-Atlas 
mountains. The natural vegetation in the Souss is savanna 
dominated by the Argan (Argania spinosa), a local endemic 
tree found nowhere else, part of the area is now a UNESCO 
Biosphere reserve to protect this unique habitat. 

The satellite images used in this study were taken respec- 
tively, in 19 March 2002 and 12 April 2005 by the Landsat 
ETM+, and are defined by the coordinates (Path 203, Row 
39). 

Figure 1 shows RGB composite of the two used images. 




(a) Landsat ETM+ 2002 




(b) Landsat ETM+ 2005 

Fig. 1: RGB composite of the Landsat ETM+ images 2002 
and 2005 (Agadir, Morocoo) 




Fig. 2: Spectral response of the topics 



C. ICM classification with constraints 

After preprocessing of the two images and establishment 
of the samples, a supervised contextual ICM classification 
with constraints [4], [5] is applied to the two images to have 
the probabilities of pixels belonging to classes. The classified 
image of 2002 is presented in figure 3. 




Fig. 3: Supervised ICM classification with constraints of the 
2002 image (Agadir, Morocco) 



D. Multidates and multi-source fusion by the hybrid DSmT 
model 

Our fusion process is composed of the following steps, first 
the definition of the framework 0, then, the estimation of the 
mass functions of each focal element by the model of Appriou, 
finally, the application of the hybrid DSmT model. 




1) Hyper power set : Taking in consideration the 

prior knowledge of the study area, we have identified 5 
classes constituting the framework 0 ( section 4.4), which 
are: Argan ( A ), Built /Oued (BO), Vegetation ( V ), 
Greenhouses ( Gh ) and Bareground (Bg). 

So, 0 is defined as follows: 0 = {A, BO , Gh , V, Bg} 

Exploiting information of the study area and also those 
obtained by ICM classifications with constraints, some ele- 
ments of the hyper power set D G seem not being adjacents 
and exclusives. To realize a better adapted study to the real 
situations, some exclusivity constraints will be taken (hybrid 
DSmT model), for example A D Gh = <p, which reduces the 
number of focal elements of the D e . 

2) Choice of threshold: The decision is made for the simple 
classes and the classes of intersection by using our rule 
decision defined previously ( section 2.5), the threshold should 
be determined in advance by experimentation and analysis 
of total uncertainty distribution after standardization, wich is 
presented in figure 4. 




Fig. 4: The standardized global uncertainty distribution of the 
fusion map 

We have tested our decision rule with various values of 
threshold, the following table I presents the occupancy rates 
of the simple classes (stable zones) and composed classes 
(change zones) according to the threshold. 



TABLE I: Occupancy rates of the simple classes and the 
composed classes according to the value of the threshold 



value of the threshold 


(%)of stable zones 


(%)of change zones 


1.0e-26 


0% 


100% 


1.0e-019 


9.64% 


90.36% 


1.0e-016 


32.12% 


67.88% 


1.0e-014 


54.408% 


45.592% 


1.0e-012 


76.86% 


23.32% 


1.0e-08 


99.99 


%0.01% 



The choice of threshold depends to the good results of 
detection with the use of the changed/ unchanged samples 
between the two dates (2002 and 2005), for this we have taken 
a suitable threshold that equal to l.Oe -14 . 

The results obtained with the various values of the threshold 
are as follows: a map contains 100% of the composed classes 



for the threshold of l.Oe— 26, 0.01% for the threshold of l.Oe— 
8 and 45.592% for the threshold of l.Oe — 14. The first result 
does not match with the ground truths because several areas of 
the region are unchanged during the considered period. The 
second result is far from the dynamic reality of the region 
that has known a great change. The last result is coherent and 
appears to be close to reality basing on the test samples of 
the stable zones and the zones of changes between the two 
images, while for thresholds which exceed the surrounding of 
(l.Oe — 14) there is a degradation of the change detection 
3) Validation of the results: The fusion map obtained with 
an adequate threshold =( l.Oe- 14) is presented in figure 5. 




Legend 
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Fig. 5: Fusion map obtained with threshold of (l.Oe- 14) 



From the fusion map, obtained with an adequate threshold 
we obtain the table II presenting the occupancy rate of the 
classes. 



TABLE II: Occupancy rate of the classes 



Classes 


(%) 


A 


29.04% 


BO 


4.869% 


Gh 


3.18% 


V 


15.92% 


Bg 


1.399% 


An BO 


0.649% 


Any 


26.40% 


AnBg 


1.87% 


BOnGh 


1.549% 


bo n v 


7.25% 


BO n Bg 


1.08% 


Ghnv 


1.98% 


Gh n Bg 


0.5% 


VHBg 


3.30‘% 



The table II illustrates the occupancy rates of the stable 
zones (simple classes) which reaches the rate of 54.408% and 
that of the zones of change (composed classes) which reaches 
45.592%. From the table II, we note that the Argan class (A) 
and Vegetation class (V) have known a great change compared 
to the other classes, Indeed, we found for the composed class 
(AnV) a rate of change of 26.40%. 

The figure 6 presents the stable zones map obtained from 
the fusion map. 




From the stable zones map, we note that the zone of Oued 
Souss (BO), the zone of the urban areas along the road Agadir- 
Taroudant (BO), the international airport of Almassira (BO), 
the parcels of the Ouled teima region (V), the greenhouses 
(Gh) of the Biogra region and also some surface of argan 
(A) are all well detected as stable zones and are assigned to 
simple classes. This attribution is well justified because some 
zones are built and are unchangeable by nature. For the change 
zones, figure 7 illustrate those zones obtained from the fusion 
map after post-processing. 




Fig. 7: Post-processed map of change zones obtained from 
fusion map 

In the post-processed fusion map, areas that have known 
major changes are the Argan class wich becomes bare ground 
or vegetation, which is well explained because of the defor- 
estation. As result we found a rate of 18.95 of the Argan class 
changed to vegetation in 2005. 

a) validation per spectral signature of the results: To 
evaluate our proposed method, we have chosen to use spectral 
signatures of classes (themes) for the two images. For pixels 
of change areas , we have a variation in the spectral signature 
of the image pixels between 2002 and 2005. Conversely, the 
pixels of the stable areas have shown a stability of spectral 
signature between the two dates. 



For this spectral evaluation, we have used the zoom window 
of the ENVM.O software for the images ETM+ 2002, ETM+ 
2005 and the fusion map to link them for locating the change 
zones in the fusion map (figure ??), Then we have compared 
the spectral signatures of the same area for the two images 
by taking in account the distributions of spectral signatures of 
different themes. 

The following figures 8 present an example of validation of 
the change results. 




(a) ETM+ 2002 (b) ETM+ 2005 (c) Fusion map 




(d) Spectral signature of V (ETM+ 
2002) 




(e) Spectral signature of A (ETM+ 
2005) 



Fig. 8: Spectral signature of the Argan class (A)(ETM+ 2002) 
and of the Vegetation class (V) (ETM+ 2005) 



From the figure 8, we note that the extract area has known 
a change of Argan theme (Class A) to vegetation theme 
(class V), what is shown by the change of the pixels spectral 
signature of the extract area which had in 2002 a spectral 
signature of the Argan (Class A), and became in 2005 that of 
the vegetation (Class V). 

The figure 9 shows the post-processed map of the fusion. 



Legend 
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Fig. 9: Post-processed fusion map 




The table III illustrates the detection of the change by fusion 

TABLE III: Rate of change obtained by fusion between 
LANDSAT ETM+ 2002 and ETM+ 2005 



Classe 


Number fo pixels 


Occupancy rate(%) 


A 


104,542 


29.0394 


BO 


17,538 


4.8717 


Gh 


11,458 


3.1828 


V 


57,320 


15.9222 


Bg 


5,048 


1.4022 


A BO 


797 


0.2214 


A^V 


68,202 


18.9450 


A- Bg 


4,295 


1.1931 


BO - Gh 


4,650 


1.2917 


BO 


25,587 


7.1075 


BO — Bg 


3,227 


0.8964 


Gh V 


3,112 


0.8644 


Gh — Bg 


471 


0.1308 


V^Bg 


2,662 


0.7394 


V - A 


26,889 


7.4692 


Bg^ A 


2,439 


0.6775 


Gh -► BO 


915 


0.2542 


V - BO 


502 


0.1394 


Bg — BO 


647 


0.1797 


V - Gh 


4,005 


1.1125 


Bg Gh 


1,324 


0.3678 


Bg^V 


12,831 


3.5642 



with :Argan(A), Built/Oued (BO), Greenhouses (Gh), Vegeta- 
tion (V), Bare ground (Bg). 

V. Conclusion 

In this paper we have proposed a new method for land 
cover changes detection. In first, we have included the spatial 
information in the process of fusion/classification by using 
a hybrid DSmT model with the introduction of contextual 
information using ICM classification with constraints, the 
use jointly of DSmT and ICM improves performance of the 
changes detection in terms of accuracy and exactitude. Sec- 
ondly, we have proposed a new decision rule that has shown 
its performance and allowed us to overcome the limitations 
of rules decision based on the maximum generalized belief 
functions which are increasing and unsuited to the decision of 
the union and the intersection of elements. 

The application of this method for the catography of the land 
cover changes is promising, however, the determination of the 
adequate decision rule (currently we work on a dynamic fusion 
method) and the addition of a temporal constraint to allow 
the processing of different dates are still significant issues for 
further investigation. 
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